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ABSTRACT 

Velocity shifts and differential broadening of radio recombination lines are 
used to estimate the densities and velocities of the ionized gas in several hyper- 
compact and ultracompact HII regions. These small HII regions are thought to 
be at their earliest evolutionary phase and associated with the youngest massive 
stars. The observations suggest that these HII regions are characterized by high 
densities, supersonic flows and steep density gradients, consistent with accretion 
and outflows that would be associated with the formation of massive stars. 

Subject headings: stars: early type — stars: formation — ISM: HII regions 
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) that are characteristically different from many of the few 
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thousand ultracompact HII (UCHII) regions scattered throughou t the galaxy (IWood fe Churchwell 



19891 : iKurtz et al.1 ll994J : IWalsh et all Il998l : iGiveon et al.l l2005al lbh. First, these very small 



HII regions have electron densities that are higher (n P > 10 5 cm 3 ) than those of UCHII 
regions (n e ~ 10 4 cm" 3 ) jKurtzl I2OO0I ; IKurtz fe Francol I2OO2I ; IPratap et al~1ll992[ ). It is un- 



clear how high the densities might be, because the density cannot be reliably determined 
from the radio continuum if the gas is optically thick. Second, when observed at centimeter 
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2006). Not all the small HII regions have 



this property, but those that do are sometimes referred to as broad recombination line objects. 
The breadth of the RRLs could be caused by a combination of electron impact (pressure) 
broadening associated with high gas density plus broadening caused by spatially unresolved 
gas motions. However, the relative contributions of broadening due to gas pressure and gas 
motions are not known. Third, these HII regions often have continuum spectral energy dis- 



tributions (SEDs) that are roughly 
the centimeter to the sub millimeter 



inear with frequency over decades in wavelength from 



Hofn er et al.lll996l : iFranco et al.ll2000l : iTesti et al.ll2000 



Beuther et al.ll2004j ; llgnace fe Churchwell! 120041 ) . The underlying cause of this scaling is not 
known although spatially unresolved density variations are a possibility. 

These small and dense HII regions are sometimes referred to as hypercompact HII 
(HCHII) regions. While the term hypercompact implies a size smaller than ultracompact 
(UC), the three properties of high electron densities, broad recombination lines, and linear 
SEDs, are also seen in some UCHII regions (~ 0.1 pc). The common characteristic between 
these unusual UCHII regions and the HCHII regions might be that they are all young and 
at the earliest evolutionary stage, but the unusual UCHII regions might contain more or 
brighter stars than their smaller HCHII counterparts. The HCHII and unusual UCHII 
regions represent some of our best sources of information on the processes in the formation 
of massive stars, and this motivates our understanding of their special properties. 

In this paper we report on multi-frequency RRL observations that address the origin 
of the three properties of the small HII regions and the processes of massive star formation. 
Because broadening due to pressure is dependent on frequency whereas broadening due to 
gas motions is independent of frequency, the widths of RRLs observed at more than one 
frequency can be used to separate the contributions of each broadening mechanism. Since 
the pressure broadening is proportional to the density, this separation allows a determination 
of the electron density. Similarly, since the broadening due to gas motions is proportional to 
the range of velocities, we can also measure the magnitude of the gas velocities. 
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2. Observations 

We selected a few well-studied, very bright, very small HII regions that previous obser- 
vations suggest have one or more of the three properties of the small HII regions. For each 
source we have one observation of the H30a line made with the Submillimeter Array (SMA) 
and one or more observations of centimeter wavelength RRLs made with the National Radio 
Astronomy Observatory Very Large Array (VLA)0 

The SMA observations of the H30a line (231.900959 GHz) in the HII regions G10.6- 
0.04, NGC7538-IRS1, W51e2, and G28.20-0.04 were made in 2005 September with a spectral 
resolution of 0.83 km s -1 , a bandwidth of 2 GHz, and an angular resolution of 1.0". Obser- 
vations of G45.07+0.14 were made in 2005 October with an angular resolution of 0.4". The 
noise level obtained in each observation was 60 mJy beam -1 channel -1 . Analysis of the data 
indicated that some of the baselines were not well-defined, hence the positions and flux den- 
sities of the sources are uncertain. The standard deviation of the flux of 3C454.4, measured 
once for each observation, is 39% of the average measured value (tabled]). If we exclude 
the low measurement of 3C454.4 corresponding to the G45.07+0.14 data, then the standard 
deviation falls to 14% of the average. Our positions from the SMA agree with those of the 
VLA to better than the angular resolution. The HII regions are unresolved by the SMA 
observations, and hence the widths of their spectral lines are unaffected by baseline errors. 
The SMA data were processed in the SMA data reduction package MIR, and in MIRIAD. 

Most of the VLA observations were made in 2003 and 2005. We reprocessed observations 



of the H76a line in G45.07+0.14 (IGaray et al.lll986l ) from the V LA ar c hives. The observa- 



tions of the H66a line in G10.6-0.4 were previously reported in iKeto J2002ah . The other 
observations have not previously been published. Details of the observations are presented 
in table [2j The VLA data were processed in AIPS; calibration data are given in table [3j The 
recombination line spectra are shown in figures 1 to 5, and the line widths and velocities 
based are given in table HI 



3. Gas velocities, pressure broadening, and electron densities 

Radio recombination lines in HII regions are broadened by three contributions: thermal 
broadening, dynamical broadening due to spatially unresolved motions in the HII region 
including both ordered flows and turbulence, and pressure broadening due to high electron 



1 The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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densities. We seek to estimate the separate contributions of each of the three broadening 
components by comparing the widths of RRL observed at different frequencies. This is 
possible because the width due to thermal and dynamical broadening, is independent of 
frequency while the width due to pressure broadening decreases with frequency as Av ~ z/~ 4 
(IBrocklehurst fc Seatonlll972l ; lGriemlll974l ). Since the dynamical broadening is proportional 
to the gas velocity while the pressure broadening is proportional to gas density, from a pair 
of RRL observations we can estimate the velocities and densities in the ionized gas. 

We determine the contributions of each of the three broadening mechanisms as follows. 
We assume that the dynamical broadening, Ai/p, and the thermal broadening, Aut, combine 
in quadrature, 

Au G = 



Avl 



Av\. 



1 



This width due to gas motions, both thermal and spatially unresolved ordered flows and 
turbulence, combines with the Lorentzian width, Avl, due to pressure broadening to produce 
a Voigt profile with width, Av v . We approximate the width of the Voigt pr ofile by an 



2002 . equation 



algebraic expression (lOlivero Sz Longbothuml 119771 : I Gordon fc Sorochenkol 
2.72), 

Av v = 0.5343 Au L + [Au G + (0.4657Az/ L ) 2 ] 1/2 (2) 

where Avy is the full-width at half-maximum (FWHM) of the Voigt profile, Au^ the FWHM 
of the Lorentzian, and Ave the FWHM of the Gaussian. Because of the rapid decrease in 
pressure broadening with frequency, the width of an RRL at a high enough frequency must 
be due solely to thermal broadening and spatially unresolved gas motions. Therefore we 
assume that this broadening, Ave in equation [2l is given by the observed width of our 
highest frequency RRL, H30a. If we compare the width of the H30a line to that of a 
lower frequency line, any increased width observed in the lower frequency line must be due 
to pressure broadening because the width Az/ G is independent of frequency We can then 
calculate the width due to pressure broadening, Avl, by solving the quadratic equation [21 

The range of gas velocities is estimated from equation ([I]). Assuming an ionized gas 
temperature of 8000 K, the dynamical widths are 15.3, 8.3, 27.1, 18.7 and 53.9 kms" 1 for 
G10.6-0.4, G28.20-0.04N, G45.07+0.14, W51e2, and NGC 7538-IRS1 respectively. 

The electron density may be estimated from the ratio of Aul and Avt- For a (AN = 
1) transitions, the ratio of the Lorentzia n width, Avt„ res ulting from pressure broadening, 
and the thermal width, Ai/?, is given by (IKeto et al.lll995h . 



Av L 
Avt 



1.2 ( _ 

vio 5 



N\ 7 
92J 



(3) 



where n e is the electron density in cm 3 , and iV is the principal quantum number. 
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As an example, we calculate the electron density in G10.6-0.4 from the H30a and H92a 
RRL observations as follows. We use the width of the H30a line, 24.5 kms -1 (table H|) to 
estimate Avq at the frequency of the H92a line. Thus Avq = 0.678 MHz. The measured 
width of the H92a line is Az/y(H92a) = 1.03 MHz. With these values and equation [21 the 
estimate of the width due to pressure broadening is Avl = 0.578 MHz or 20.9 kms -1 . The 
thermal width, assuming a temperature of 8000 K, is Aut = 1.43 MHz or 19.1 kms -1 . The 
width due to gas motions is Av^ = V24.5 2 — 19. 1 2 = 15.3 kms -1 . The electron density 
from equation [3] is 9.1 x 10 4 cm -3 . Similarly, from the observations of the other centimeter 
RRLs observed in G10.6-0.4, the electron densities are 7.8 x 10 5 cm -3 from the H66a line 
and 2.5 x 10 6 cm -3 from the H53a line. Comparing these three RRL, the electron density 
increases with the RRL frequency. If there is a density gradient in G10.6-0.4 this would 
suggest that the higher frequency lines are seeing deeper into the nebula where the electron 
densities are higher. 

3.1. Uncertainties 

We can check our approximation of negligible pressure broadening in the H30a line from 
[3j For the H30a line, this ratio is less than one as long as the density is less than 3 x 10 s 
cm -3 . Thus our assumption of negligible pressure broadening in the H30a line is valid up 
to this density. 

The lines in our sample are not consistently Gaussians or by Lorentzian. One possibility 
is that the asymmetric line profiles indicate asymmetries in the gas motions in the HII regions, 
for example, ordered flows rather than isotropic turbulence. We determine the line widths 
by fitting Gaussians because this provides a reasonable estimate if the observed profiles differ 
from a Gaussian most strongly in the wings and are similar in the core where the width is 
measured. A comparison of the line widths derived by fitting Gaussians to Voigt profiles 
provides an informal estimate of the possible error introduced by this method. If the Voigt 
profile is composed of a Gaussian and a Lorentzian of equal FWHM, then a Gaussian fit 
overestimates the true FWHM of the Voigt by 10%. 

We find that the procedure of using Gaussian line widths in equation [2] to estimate the 
pressure broadening is a more reliable method than some alternatives. For example, one 
could imagine fitting a Voigt profile directly to an observed low frequency line to determine 
simultaneously both the Gaussian and Lorentzian components. However, we find this unreli- 
able because the fitting procedure assigns the line width to either the Gaussian or Lorentzian 
components primarily on the shape of the line wings where the profiles differ most strongly, 
but the signal-to-noise ratio is the weakest. Alternatively, one might estimate the FWHM 
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directly from the 3 channels with the peak and half-power emission. However, because this 
estimate relies on single channel measurements it is less reliable than a line-fitting procedure 
that uses measurements at all the channels across the line profile. 

Several lines of evidence, discussed later in the paper suggest that there are density 
variations within the HII regions. Because all the HII regions in our sample are unresolved, 
the measured electron density represents the average within the HII region. The size of 
the HII region within the beam (filling factor) does not affect this estimate of the average 
density because the electron densities are determined from the line widths rather than the 
emission intensity. However, the unknown density structure of the HII region creates some 
uncertainty. If there are significant variations in optical depth with frequency then it is 
possible that lines at different frequencies could be generated in different locations within 
the HII region. 

Aside from the systematic effects, the uncertainty in the average electron density is given 
by the standard propagation of errors, assuming the errors in the widths are independent. 
Uncertainties in the electron densities are listed in table and for the line center velocities 
and widths in table HI 



T 

t v = 8.235 x 1(T 2 ( 



(4) 



3.2. Electron densities from the radio continuum 

At the high densities in our HII regions, the centimeter continuum emission is opti- 
cally thick to free-free em ission. This can be det ermined from the formula (eqn. A. lb, 
Mezger fc Henderson! 1 196 71 ) or (eqn. 10. iKetol 120031 ) . 

1.35, v y2.l, n 2 L 

KJ VGHzJ Vpc cm- 

where T e is the electron temperature and n 2 e L is the electron density squared times the path 
length. For a characteristic size of 0.01 pc, an HII region has an optical depth of unity at 
22 GHz (1.3 cm) if the elec tron density is 4 x 10 5 cm" 3 . This means that the formulae for 
optically thin emission (e.g. iMezger fc Henderson! 119671 ) are not useful for the HII regions in 
our sample, and probably not for small, high-density HII regions in general. 

If the gas is optically thick to free-free emission, and the electron density is estimated 
from the radio continuum by incorrectly assuming optically thin emission, then the derived 
electron density is just that density required to make the optical depth unity. This can be 
understood by a simple example. The radio free-free emission through a slab of uniform gas 
density and temperature is, 



S v = 2kT e (v 2 /c 2 )[l - e-^KJy/Sr) 



(5) 
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The optically thin approximation amounts to replacing [1 — exp(— t v )\ by r v so that, 

S v « 2kT e (u 2 /c 2 )r u (6) 

whereas, in the optically thick limit, 

S„ « 2kT e {u 2 /c 2 ) (7) 

Thus if the HII region is optically thick (equation [7j) but r v is derived assuming equation [BJ 
then the derived value of the optical depth is 1, and the electron density derived from equation 
IU is the density required to make r u = 1. Thus at 22 GHz (1.3 cm), the electron density 
derived for an HII region with a size of 0.01 pc is never be greater than about 10 5 cm -3 if 
the optically thin formulae are used - regardless of the true density. This underestimate also 
affects quantities derived from the electron density such as the number of ionizing photons 
required to balance recombinations which in turn is often used to derive the spectral types 
of the exciting stars. 

Formulas for spherical HII regions that have a density grad ient such that the central 



part of the HII region is optically thick are derived in iKetol (120031 ). However, these formulae 



require a specification of the size of the optically thick region. 



3.3. The correlation of line widths, velocities, and frequencies 



As shown in table HI the line widths decrease with increasing frequency in all cases. 
This suggests the increasing importance of pressure broadening at lower frequencies. For 
sources with more than one centimeter RRL observed, the electron densities increase with 
the increasing frequency of the centimeter line used in the analysis. The one exception is the 
electron density calculated from the H66a line in NGC7538-IRS1. This line is very much 
wider and out of character with the other measurements in our sample. Aside from this 
one exception, the correlation of density with observing frequency indicates that the higher 
frequency lines are associated with higher density gas. The line center velocities also increase 
(redshift) with frequency, except for the lines of NGC7538-IRS1. 



The combination of the two correlations, line center velocity wi t 



(1987 


) and 


Keto et al. 


(1995) 



i frequency and line 



ing flow within a density gradient (IGuilloteau et al.l 119831) . A full explanation is given in 



Keto et al. ( 1995 ) (see also Brocklehurst fc Seaton ( 1972l ))and can be summarized as follows: 



The peak intensity of RRLs is proportional to the density squared and inversely proportional 
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to the line width. Thus the peak intensity of optically thin high frequency lines with negligi- 
ble pressure broadening is proportional to the density squared. The peak intensity of lower 
frequency lines is proportional to the first power of the density since the pressure broadened 
line width is itself proportional to density. Thus in spatially unresolved observations of gas 
with a density gradient, the higher frequency RRLs are dominated by emission from the 
higher density gas while the lower frequency RRLs include contributions from both high and 
low density gas. An observed difference in velocity between the high and low frequency lines 
thus indicates a difference in velocity between the high and low density gas. 

The correlation of velocities and width with frequency therefore indicates that in these 
HII regions there is both a flow and a density gradient. The two are naturally related by the 
conservation of mass and the geometry of the flow. In general, if the flow is either outward 
and diverging or inward and converging then conservation of mass requires that there be 
a density gradient except in the unique circumstance that the variation of the flow speed 
exactly cancels the geometric divergence or convergence. 



4. Continuum Spectral Indices 



We determined the spectral energy distrubutions (SEDs) of the HII regions in our sam- 
ple using our own millimeter and centimeter observations and data collected from the litera- 
ture. We selected only those data with angular resolutions close to 1" specifically excluding 
single dish observat ions. Most o f the data fo r NGC 7538-IRS1 and G45.07+0.14 are from 
Pratap et al.l (119921 . table 1) and iGaray et al.l (119861 . table 3), respectively. The rest of the 
observations are listed in table Most of these data do not include estimated uncertainties. 
However, the calibration of radio data is fairly standard. Many of the reported observations 
were made at the VLA that typically calibrated flux densities to better than 20% at 1.3 cm 
at the dates of the observations and much better at longer wavelengths. As noted in §21 we 
estimated our SMA observations to be calibrated to better than 39%. 

The spectral energy distributions are shown in figures 6 through 10. Along with each 
observed SED, shown as crosses, the figures show two model SEDs as solid lines. One model 
SED is based on a single HII region with a power- law density gradient. The continuum 
flux density for this model is calculated as in iKetol (120031 ). The radio free-free emission of a 
spherical HII region seen as an unresolved source is, 



S u = A7ckT e u 2 /c 2 / 6[1 
Jo 



- T »]d9 



Here, 6 is an angular coordinate on the plane of the sky, 6 = ■sjx 2 + y 2 / D with coordinate 
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axes x and y are centered on the HII region, 8q = Rq/D where Rq is the radius of the HII 
re gion, and D the distance to the source. An approximation to the free-free optical depth 



is. 



Mezger fc Henderson! (119671 . equation A. la), 



r v = 8.235 x 10" 



-1.35 



V 



-2.1 



EM 



J \GHzs \pc cm 
where T e is the electron temperature, and EM is the emission measure, n 2 e L. 



(9) 



The model parameters for all sources are listed in table [71 The other model is an SED 
based on constant density HII regions. The figures for NGC7538 - IRS1 and W51e2 show the 



SEDs from multi-component models suggested by iPratap et al.l (119921 ) and iRudolph et al. 



(Il990l ). T hese are both based on constant density HII regions plus dust. The dust emission 
enters as (IPratap et al.lll992l . equation 1), 



S„(Jy) = 52.36M 



0.2 \^ 



exp 



14.4 



1 



-i 



^kpc 



(10) 



where S v is the dust emission in Jy, M is the mass of the dusty gas in M Q , (3 is the spectral 
index of the dust emissivity, is the dust temperature in Kelvin, and D is the distance in 
kpc. 



The multi-component model for NGC7538-IRS1 (IPratap et al.lll992l ) nicely illustrates 
how a density gradient within an HII region can stretch the transition between optically 
thick and thin frequencies resulting in an SED with an intermediate spectral index. In this 
model there are two concentric spherical regions each of uniform but different densities. The 
smaller, higher density component produces an emission curve that is shifted upward and 
to the right of the curve for the larger, lower density component (figure 6). The emission 
from the two components combines to produce an SED of intermediate slope. An HII region 
with a continuous density gradient can be considered as the log ical extension of this two 
component model to many nested smaller and denser HII regions. iLugo et al.l ( 2004(@f ) also 
model the NGC 7538-IRS1 HII region with a density gradient that they suggest is due to a 
wind off a photo-evaporating disk. 

The SED modeling combined with the recombination line measurements shows that dust 
is not a significant contributor to t he cont i nuum e mission measured at high angular reso- 



lution from these sources. (See also iKurta (120051 ) . iDePree et al.l (119981 ) . and iDePree et al. 
( 120001 ).) Consider W51e2 as an example. As seen in figure 7, the HII region can be mod 



eled with u niform, low- density ionized gas, with the 147 and 230 GHz emission arising 
from dust (IRudolph et al.l Il990l ). Figure 7 shows that with this approach, the free- free 
continuum emission at 230 GHz would be about a factor of 10 below the observed emis- 
sion. However, the line-to-continuum ratio should be approximately one at this frequency 
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( IGordon fc Sorochenkoll2002l ). Hence, if the free- free continuum were this low, the H30a re- 
combination line would also have a much lower flux density — so low that the line would not 
be detected in our observations. Stated another way, the recombination line emission cannot 
be 10 times greater than the free-free continuum emission at this frequency. Therefore, most 
of the continuum emission observed by the interferometer at 231.901 GHz must be free-free, 
hence the electron density must be higher, and a density gradient is required. 

The high frequency emission observed in large beams (20 — 30") ar ound these bright 
HII regions certainly arises from dust. Again take W51e2 as an example. iJaffe et al.l (119841 ) 
report a flux densit y of 1200 Jy beam -1 a t 750 GHz and 40" angular resolution. We assume, 
as in the model of iRudolph et al.l (119901 ) . that this emission is from cool dust in the large- 
scale molecular cloud surrounding the HII region. Because th e dust is cool (50 - 100 K), a 



large mass of dust is required (4500 M , IRudolph et al.l (Il990h ) to produce the observed 750 



GHz flux density. While the emission from this dust is also sufficient to produce the 90 and 
230 GHz flux density that is observed by interferometers, the cool dust is not the source of 
this emission because it is not possible that this large mass of cool dusty gas is within the 
spatial scale defined by the arc second angular resolution of the interferometers. However, 
the ionized gas is sufficiently hot (10 4 K) that only a small mass (0.01 M ), consistent with 
the mass of the HII regions, is required to produce the 90 - 230 GHz flux density. Similar 
considerations apply to the multi-compone nt models of other H II regions, for example, the 
model discussed earlier for NGC7538-IRS1 (IPratap et al.lll992l ). 

On a spatial scale similar to that of HCHII and UCHII regions, continuum emission 
from dust has bee n observed by the S MA in the nearby (1.7 kpc) massive star forming re- 
gion NGC6334 by iHunter et al.l (120061 ). The dusty clumps of arcsecond scale in NGC6334 
have continuum emission on the order of 1 Jy. Scaled to the distance of W51, the dust in 
such clumps would be a minor contributor (60 mJy) to the continuum emission. Thus, high 
angular resolution observations of dusty clumps in NGC6334 are consistent with the inter- 
pretation of negligible dust contribution to the emission observed at high angular resolution 
in more distant bright HII regions. 



5. Discussion 



The observations suggest that pressure broadening of the centimeter RRLs is impor- 
tant in small HII regions. This contradicts the interpretation of some previous studies. In 
particular, on the basis of e lectron densit i es calc u lated by assuming c onsta nt gas density 
and o ptically thin emission, iGaume et al.l (Il995al ). iDePree et al.l (119971 ). and iDePree et al. 
(120051 ) have suggested that the electron densities in NGC7538-IRS1, W49, and SgrB2 are 
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too low to cause significant pressure broadening in the H66a line. However, the RRL ob- 
servations suggest that the electron d ensity in NGC7538-IR S1 is two orders of magnitude 



greater than the 10 5 cm 3 derived by iGaume et al. 



(1995a) and sufficient to cause signifi- 



cant pressure broadening. iJaffe fc Martm-Pintadd (119991 ) report observations of millimeter 
RRL of several sources, including G10.6-0.4 and NGC7538-IRS1. They concluded that the 
RRL widths are not due to pressure broadening because the observed widths of different 
recombination lines do not scale as the fourth power of the frequency. Their conclusion 
was based on a comparison of observations with different angular resolutions. It has been 
demonstrated that low ang ular resolution o b servations may suffer confusion from more- 
extended, lower-density gas ( Simpson 1973a b; Gulvaey fc Sorochenko 1973 ; Smirnov et al. 



984; Gordon &; Sorochenko 



2002). We co ncur with IGaume et al.l (I1995af ): 



DePree et al. 



(119971 . 120051 ) ; I Jaffe fc Martin-Pintado! (119991 ) that broadening due to gas motions is a signif- 



icant contribution to the observed line widths, but we find that pressure broadening is also 
significant, particularly for the lower frequency lines. 

The intermediate SEDs of small HII regions imply a density gradient as would arise from 
a divergence or convergenc e in th e flow of ionized g a s in a n HII region. Oth e r inte rpretations 
are possible. Franco et al.l (120001 ) . iKurtz fc Franco! (120021 ) and iKim fc Kool (120021 ) suggested 
that the various densities are due to the hierarchical structure of the molecular cloud prior 
to its ionization by a massive s tar. The intermediate SEDs can also be created by density 
variations other than gradients. Ilgnace fc Churchwelll (120041 ) suggested that an HII region is 
actually an ensemble of differently-sized clumps, each of uniform but different density. One 
problem with these models of pre-existing structure or clumps is that their survival time is 
approximately their sound crossing time. This timescale seems impossibly short if the high 
density clumps are substructures within an HCHII (t c < 1000 yrs for R < 0.01 pc). 

The interpretation of density gradients and resulting high electron densities in small 
HII regions also contradicts the interpretation of significant dust emis s ion in a number o f 
other studies For e xample. iRudolph et al.l ffl99ph : IPratap et all Jl992h : iTesti et all (boooh : 
Schilke et al.Ul99oh : lBeuther et all (hooibood ) model the small HII regions, W51e2, NGC753£ 
IRS1, G9.62+0.19-F, and Orion-KL-SMAl respectively, with low a nd constant densit y gas 
and assume that dust makes up the emission at high frequencies. iTesti et al.l (120001 ) also 
provide an alternative model with a stellar wind rather than dust to supply the high fre- 
quency emission. These models with low and constant density have difficulty in explaining 
the broad RRL line widths. If the density is low, then the widths must be due to dynamical 
broadening created by a supersonic flow. However, a constant density is not generally con- 
sistent with supersonic flows. In general, a strong density gradient is required to create the 
strong pressure gradient needed to drive the flow. 
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6. Conclusions 

1. The radio recombination line widths increase with decreasing frequency in all cases. 
This indicates the greater importance of pressure broadening at lower frequencies. 

2. The electron densities calculated by comparison of centimeter recombination lines 
with the H30o; millimeter line increase with increasing frequency of the centimeter line. This 
suggests that the higher frequency lines are generated in denser gas and suggest that density 
gradients exist in small HII regions. 

3. Supersonic line widths are observed in the H30a RRL. This line is at a frequency 
high enough that pressure broadening is negligible provided that the density of the emitting 
gas is less than 10 8 cm -3 . The large line widths suggest the presence of supersonic motions 
within the HII regions. 

3. The electron densities in HCHII and UCHII regions are generally high enough that the 
gas is optically thick at centimeter wavelengths. This means that electron densities cannot 
be calculated from radio continuum observations using formulae that assume optically thin 
emission. Instead formulae appropriate for partially optically thick emission must be used. 

4. At high frequencies (> 100 GHz) and high angular resolution (~ 1") the observed 
continuum emission is mostly free- free emission and not from dust. Comparisons based on 
high frequency observations with different beam sizes can lead to erroneous conclusions. High 
frequency emission at low angular resolution (> 10") is probably from dust, but associated 
with the large-scale, overlying molecular cloud. 

5. Kurtz acknowledges support from UNAM, DGAPA project IN106107. 
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Table 1. Calibration of Millimeter Observations 







Flux 




Bandpass 




Gain 




Source 


Band 


Calibrator 


Flux 


Calibrator 


Flux 


Calibrator 


Flux 








Jy 




Jy 




Jy 


G10.60 -0.4 


SMA 


Ceres 


0.7 


3C454.3 


16.8 


PKS 1911-201 


1.4 


G28.20 - 0.04 N 


SMA 






3C454.3 


14.3 


PKS 1911-201 


1.8 


G45.07 + 00.14 


SMA 






3C454.3 


4.6 


PKS 2025+337 


0.5 


W51e2 


SMA 


Ceres 


0.7 


3C454.3 


16.0 


PKS 2025+337 


0.8 


NGC7538-IRS1 


SMA 


Ceres 


0.5 


3C454.3 


11.9 


PKS 2202+422 


3.0 



Note. — Flux densities derived assuming a flux density for Uranus of 36.6 Jy. 



Table 2. VLA Observations 



Source 


RA 


Dec 


Line 


Frequency 


Spectral 


Band 


Angular 


Cont. a 


Date 












Resolution 


width 


Resolution 


rms 


Observed 










GHz 


kms 


Chan. 


arc second 


mjy 




G10.6 -0.4 


18 10 28.7 


-19 55 49 


H92a 


8.309383 


1.76 


64 


1.5 X 1.0 


7 


1989-03-17 








H66a b 


22.36417 


2.62 


64 


1.7 x 1.0 


9 


1988-05-23 








H53a 


42.95197 


1.36 


64 


1.5 x 1.3 


3.0 


2003-02-03 


G28.20 - 0.04 N 


18 42 58.17 


-04 13 57.0 


H53a 


42.95197 


1.36 


64 


0.6 x 0.4 


1.2 


2005-06-24 








H53a 


42.95197 


5.45 


32 


0.7 x 0.5 


3.0 


2005-06-28 


G45.07+0.14 N 


19 13 22.069 


10 50 52.5 


H76a c 


14.69000 


15.95 


16 


0.6 x 0.6 


10 


1984-01-07 


W51e2 


19 23 43.913 


14 30 14.7 


H66a 


22.36417 


2.62 


64 


0.3 x 0.2 


1.6 


2003-02-01 








H53a 


42.95197 


1.36 


64 


0.5 x 0.4 


0.9 


2003-02-01 


NGC7538-IRS1 


23 13 45.37 


61 28 10.4 


H53« 


42.95197 


5.45 


31 


0.2 x 0.2 


0.9 


2005-04-17 



a Bandwidth = channels X spectral resolution 

Data from the VLA archives. The original observation is reported in iKetd d2002af) . 
c Data from the VLA archives. The original observation is reported in lGarav et al. | lll986h . 



Table 3. Calibration of Centimeter Observations 







Flux 




Bandpass 




Gain 




Source 


Band 


Calibrator 


Flux 


Calibrator 


Flux 


Calibrator 


Flux 








Jy 




Jy 




Jy 


G28.20 - 0.04 N 


Q 


3C286 


1.5 


3C273 


24.2 


PKS 1851+005 


1.9 


G45.07 + 00.14 N 


u 


3C286 


3.5 


3C454.3 


7.5 


PKS 1923+210 


2.1 


W51e2 


Q 


3C286 


1.5 


3C273 


10.1 


PKS 1923+210 


2.6 


W51e2 


K 


3C286 


2.5 


3C273 


19.5 


PKS 1923+210 


3.1 


NGC7538-IRS1 


Q 


3C48 


0.5 


3C454.3 


15.6 


PKS 2250+558 


0.8 
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Table 4. Recombination Line data 



Source 


Continuum 
Peak 
Jy/beam 


Integrated 
Intensity 

Jy 


Quantum 
Number 


Line 
Peak 
m Jy/beam 


lo- 
rn Jy/beam 


Line 
Velocity 8, 
kms -1 


la 
kms 


Line 
Width b 
kms 


la 
kms 


G10.6-0.4 


0.017 


3.14 


92 


8.1 


0.2 


0.6 


0.5 


37.5 


1.2 




0.04 c 


4.57 


66 


188.0 


5.0 


2.1 


0.005 


35.1 


0.005 




2.1 


6.0 


53 


219.0 


4.9 


2.2 


0.005 


31.8 


0.01 




1.17 


2.41 


30 


1040.0 


28.5 


3.4 


0.3 


24.5 


0.8 


G28.20-0.04 N 


0.4 


1.1 


53 


99.0 


3.6 


88.7 


0.6 


33.4 


1.5 




0.56 


0.72 


30 


930.0 


22.0 


92.5 


0.2 


20.9 


0.6 


G45.07+0.14 


0.37 


0.77 


76 


21.0 


2.5 


45.7 


3.4 


56.9 


8.1 




0.08 


0.17 


30 


92.0 


10.2 


49.4 


1.8 


33.2 


4.2 


W51e2 


0.15 


0.30 


66 


14.2 


0.4 


53.8 


0.7 


50.9 


1.6 




0.17 


2.3 


53 


42.0 


1.1 


59.5 


0.3 


32.5 


1.2 




1.67 


3.66 


30 


870.0 


51.0 


59.7 


0.8 


26.8 


1.9 


NGC7538-IRS1 


0.2 


0.7 


53 


23.3 


1.1 


-41.8 


1.0 


61.0 


4.1 




2.56 


2.79 


30 


4450.0 


58.0 


-60.5 


0.4 


57.3 


0.8 



a Line center velocity from Gaussian fit. 
b FWHM from Gaussian fit. 
c Keto 2002 



Table 5. Electron Densities from line widths 



Source 


Quantum Number 


Line Width 


la 


n e 


lo- 






kms -1 


kms -1 


cm~ 3 


om " 3 


G10.6-0.4 


92 


37.5 


1.2 


9.1 X 10 4 


7.3 x 10 3 




66 


35.1 a 


0.005 


7.8 x 10 5 


3.2 x 10 2 




53 


31.8 


0.01 


2.5 x 10 6 


3.2 x 10 3 


G28.20-0.04 N 


92 


74.4 b 


2.6 


2.9 x 10 5 


1.2 x 10 4 




53 


33.4 


1.5 


4.1 x 10 6 


4.3 x 10 5 


G45.07+0.14 


76 


56.9 d 


8.1 


6.1 x 10 5 


1.8 x 10 5 




66 


42.3 C 


2.3 


6.9 x 10 5 


1.6 x 10 5 


W51e2 


66 


50.9 


1.6 


1.6 x 10 6 


9.0 x 10 4 




53 


32.5 


1.2 


2.0 x 10 6 


4.0 x 10 5 


NGC7538-IRS1 


66 


180 c 




1.0 x 10 7 






53 


61.0 


4.1 


1.4 x 10 6 


1.5 x 10 6 



Note. — Electron density, n e , required to produce the observed increase in the width of 
the cm RRL over the mm RRL by pressure broadening. 

a Kcto 2002 

b Scwilo ct al 2004 

c Gaume ct al 1995a 

d Garay ct al 1986. The line width measured in our analysis of this data is larger than 
reported in the original reference (48.1 kms -1 ). 



c Garay ct al 1985 
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Table 6. Integrated Intensities from Previous Observations 



Source 


Frequency 


Flux 


Beam 


Reference 




GHz 


Jv 


Arc seconds 




G10.6-0.4 


23 


2.4 


0.1 


Sollins & Ho (2005) 




22.3 


4.1 


0.2 x 0.1 


Keto & Wood (2006) 




15 


2.8 


0.3 


Ho & Haschick (1981) 




15 


2.6 


0.3 x 0.2 


Turner & Matthews (1984) 




8.3 


3.1 


1.5 x 1.0 


this paper 




5 


1.7 


1 


Ho & Haschick (1981) 


G28. 20-0.04 N 


23.7 


0.98 


0.3 x 0.2 


Sollins et al. f 2005b) 




15 


0.54 


0.5 


Kurtz et al. ( 1994) 




8.4 


0.30 


0.9 


Kurtz et al. ( 1994) 




8.3 


0.30 


4.9 x 2.7 


Sewilo et al. f2004) 


G45.07+0.14 


345 


2.71 


3.0 x 4.0 


Su et al. (2004) 




98 


1.04 


2.6 x 2.4 


Hunter et al. (1997) 


W51e2 


147 


3.7 


2.0 


Zhane et al. (1998) 




88.2 


1.23 


6x5 


Rudolph et al. (1990) 




23.7 


0.35 


5x4 


Rudolph et al. (1990) 




14.7 


0.26 


4.1 x 3.9 


Garav et al. (1985) 




8.3 


0.07 


2.4 


Mehrineer (1994) 




5.0 


< 0.03 


4.1 


Mehringer (1994) 




5.0 


0.025 


0.5 


Rudolph et al. (1990) 



Table 7. Model HII regions for continuum SEDs 



Constant Density Density Gradient 



Source 


Distance 


radius 


density 


radius 


density 


exponent 




kpc 


pc 


cm -3 


pc 


cm -3 




G010.6 - 0.4 


6.0 


0.05 


6.0 x 10 4 


0.05 


3.5 x 10 4 


-1.4 


G28.20 - 0.04 N 


9.1 


0.03 


9.0 x 10 4 


0.043 


1.0 x 10 4 


-2.5 


G045.07 + 00.14 


9.7 


0.025 


1.6 x 10 5 


0.05 


1.5 x 10 4 


-1.7 


W51e2 a 


7.0 


0.01 


3.0 x 10 5 


0.01 


1.0 x 10 5 


-2.5 


NGC7538-IRSl b 


3.5 


0.02 


2.7 x 10 6 


0.08 


2.7 x 10 6 


-2.0 


NGC7538-IRS1 C 




0.08 


5.0 x 10 4 









a SED also includes a contribution from dust modeled as 4500 Mq of cold, dusty gas with 
a tem perature of 100 K and a dust emissivity, (3 = 1.5 from the model of iRudolph et all 
lll99(]|) . 

b Constant densit y model with 3 components, 2 HII regions plus dust from the model of 
IPratap et al. I lll992l) . The first HII component SED also includes a contribution from dust 
modeled as 53 Mq of cold, dusty gas with a temperature of 50 K and a dust emissivity, 
& = 1.0 

c Second component of constant density model. 
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Fig. 1. — H66a (dashed blue line), H53a (dotted green line), and H30a (solid red line) in 
W51e2. The flux densities of the H66a line and the H53a lines have been multiplied by 80 
and 27, respectively. 



G28.20-0.04 N 
1.2 [ — 1 1 ' 1 1 1 ' 1 1 ' ' r 




-0.2 I , , , I , , , I , , , I , , , I , , , 

40 60 80 100 120 140 

velocity (km/s) 



Fig. 2.— H53a (dashed blue line) and H30a (solid red line) in G28.20-0.04N. The flux 
density of the H53a line has been multiplied by 10. 
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NGC7538-IRS1 
5f — I — ' — ' — ' — ' — I — ' — ' — ' — ' — I — ' — 




-150 -100 -50 50 

velocity (km/s) 



Fig. 3. — H30a (solid red line) and H66a (dashed blue line) in NGC7538-IRS1. The flux 
density of the H66a has been multiplied by 190. 




Fig. 4.— H76a (dashed blue line) and H30a (solid red line) in G45.07+0.13. The flux 
density of the H76a line has been multiplied by 5. 
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Fig. 5. — H92a (long-dashed green line), H66a (short-dashed blue line), H53a (dotted black 
line), and H30a (solid red line) in G10.6-0.04. The flux densities of the H92a, H66a, and 
H53a lines have been multiplied by 115, 5.5, and 4.8 respectively. 
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NGC7538 




10.0 10.5 11.0 

log frequency (GHz) 



12.0 



Fig. 6.— SED NGC 7538 IRS1. Data are from table 1 of IPratap et all (11992f ) with the 
addition of our 230 GHz SMA observation. The blue lines show a 3 component model 



proposed by IPratap et al.l (jl992j) consisting of two HII regions plus dust (dotted line). The 
dash-dotted line (of the smaller, higher density component) is shifted upward and to the 
right of the dashed line (of the larger, lower density component. The wavy solid blue line 
shows the sum of the three components. The red solid line that crosses smoothly through the 
waves of the three component model shows the SED of a single HII region with a power-law 
density gradient. See $3] and table 7. 
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W51e2 




log frequency (GHz) 



Fig. 7. — Measured flux densities of the W5 1e2 HIT region ( crosse s). The blue dashed lines 
show a two- component model proposed by iRudolph et al.l (119901 ) consisting of a uniform 
density HII region plus dust (dotted line). The solid red line is the SED predicted for an 



HII region with a density gradient 
table 7. 



n, 



oc r L5 ) and no significant dust emission. See $3] and 



GAL 010.60-00.40 




log frequency (GHz) 



Fig. 8. — Measured flux densities of the G10.6-0.4 HII region (crosses). The blue dashed line 
is a model based on a uniform density HII region; the red line is a model based on an HII 
region with a density gradient. 
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GAL 045.066+00.138 
1.01 




log frequency (GHz) 

Fig. 9. — Continuum spectral energy distribution for G45. 07+00. 14. The crosses mark the 
observations; the blue dashed line is a model based on a uniform density HII region; the red 
line is a model based on an HII region with a density gradient. 



GAL 028.20-00.04 N 




12.0 



log frequency (GHz) 

Fig. 10. — Continuum spectral energy distribution for G28.20-0.04. The crosses mark the 
observations; the blue dashed line is a model based on a uniform density HII region; the red 
line is a model based on an HII region with a density gradient. 



